home *** CD-ROM | disk | FTP | other *** search
/ IRIX Base Documentation 2002 November / SGI IRIX Base Documentation 2002 November.iso / usr / share / catman / p_man / cat3 / SCSL / chpgst.z / chpgst
Encoding:
Text File  |  2002-10-03  |  3.8 KB  |  133 lines
  1.  
  2.  
  3.  
  4. CCCCHHHHPPPPGGGGSSSSTTTT((((3333SSSS))))                                                          CCCCHHHHPPPPGGGGSSSSTTTT((((3333SSSS))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      CHPGST - reduce a complex Hermitian-definite generalized eigenproblem to
  10.      standard form, using packed storage
  11.  
  12. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  13.      SUBROUTINE CHPGST( ITYPE, UPLO, N, AP, BP, INFO )
  14.  
  15.          CHARACTER      UPLO
  16.  
  17.          INTEGER        INFO, ITYPE, N
  18.  
  19.          COMPLEX        AP( * ), BP( * )
  20.  
  21. IIIIMMMMPPPPLLLLEEEEMMMMEEEENNNNTTTTAAAATTTTIIIIOOOONNNN
  22.      These routines are part of the SCSL Scientific Library and can be loaded
  23.      using either the -lscs or the -lscs_mp option.  The -lscs_mp option
  24.      directs the linker to use the multi-processor version of the library.
  25.  
  26.      When linking to SCSL with -lscs or -lscs_mp, the default integer size is
  27.      4 bytes (32 bits). Another version of SCSL is available in which integers
  28.      are 8 bytes (64 bits).  This version allows the user access to larger
  29.      memory sizes and helps when porting legacy Cray codes.  It can be loaded
  30.      by using the -lscs_i8 option or the -lscs_i8_mp option. A program may use
  31.      only one of the two versions; 4-byte integer and 8-byte integer library
  32.      calls cannot be mixed.
  33.  
  34. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  35.      CHPGST reduces a complex Hermitian-definite generalized eigenproblem to
  36.      standard form, using packed storage. If ITYPE = 1, the problem is A*x =
  37.      lambda*B*x,
  38.      and A is overwritten by inv(U**H)*A*inv(U) or inv(L)*A*inv(L**H)
  39.  
  40.      If ITYPE = 2 or 3, the problem is A*B*x = lambda*x or
  41.      B*A*x = lambda*x, and A is overwritten by U*A*U**H or L**H*A*L.
  42.  
  43.      B must have been previously factorized as U**H*U or L*L**H by CPPTRF.
  44.  
  45.  
  46. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  47.      ITYPE   (input) INTEGER
  48.              = 1: compute inv(U**H)*A*inv(U) or inv(L)*A*inv(L**H);
  49.              = 2 or 3: compute U*A*U**H or L**H*A*L.
  50.  
  51.      UPLO    (input) CHARACTER
  52.              = 'U':  Upper triangle of A is stored and B is factored as
  53.              U**H*U; = 'L':  Lower triangle of A is stored and B is factored
  54.              as L*L**H.
  55.  
  56.      N       (input) INTEGER
  57.              The order of the matrices A and B.  N >= 0.
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. CCCCHHHHPPPPGGGGSSSSTTTT((((3333SSSS))))                                                          CCCCHHHHPPPPGGGGSSSSTTTT((((3333SSSS))))
  71.  
  72.  
  73.  
  74.      AP      (input/output) COMPLEX array, dimension (N*(N+1)/2)
  75.              On entry, the upper or lower triangle of the Hermitian matrix A,
  76.              packed columnwise in a linear array.  The j-th column of A is
  77.              stored in the array AP as follows:  if UPLO = 'U', AP(i + (j-
  78.              1)*j/2) = A(i,j) for 1<=i<=j; if UPLO = 'L', AP(i + (j-1)*(2n-
  79.              j)/2) = A(i,j) for j<=i<=n.
  80.  
  81.              On exit, if INFO = 0, the transformed matrix, stored in the same
  82.              format as A.
  83.  
  84.      BP      (input) COMPLEX array, dimension (N*(N+1)/2)
  85.              The triangular factor from the Cholesky factorization of B,
  86.              stored in the same format as A, as returned by CPPTRF.
  87.  
  88.      INFO    (output) INTEGER
  89.              = 0:  successful exit
  90.              < 0:  if INFO = -i, the i-th argument had an illegal value
  91.  
  92. SSSSEEEEEEEE AAAALLLLSSSSOOOO
  93.      INTRO_LAPACK(3S), INTRO_SCSL(3S)
  94.  
  95.      This man page is available only online.
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.